From Big Crunch to Big Bang with AdS/CFT 
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The AdS/CFT correspondence is used to describe five-dimensional cosmology with a big crunch 
singularity in terms of super- Yang-Mills theory on IR x S 3 deformed by a potential which is un- 
bounded below. Classically, a Higgs field in the dual theory rolls to infinity in finite time. But since 
the S 3 is finite, the unstable mode spreads quantum mechanically and the singularity is resolved 
when self-adjoint boundary conditions are imposed at infinity. Asymptotic freedom of the coupling 
governing the instability gives us computational control and the quantum spreading provides a UV 
cutoff on particle creation. The bulk interpretation of our result is a quantum transition from a big 
crunch to a big bang. An intriguing consequence of the near scale-invariance of the dual theory is 
that a nearly scale-invariant spectrum of stress-energy perturbations is automatically generated in 
the boundary theory. We comment on implications for more realistic cosmologies. 



A fundamental challenge facing string and M-theory is 
that of dealing with the big bang singularity. Not only 
did the visible cosmos emanate from it, so too did the ex- 
tra dimensions whose configuration determines the laws 
of particle physics and key parameters such as the dark 
energy density. Hence it is vital both to particle physics 
and to cosmology that we resolve the initial singularity. 

It is natural to turn for guidance to the AdS/CFT 
correspondence, []J which has emerged as an extremely 
powerful tool for studying quantum gravity. According 
to AdS/CFT, a gravitational theory in an asymptotically 
Anti-de Sitter (AdS) spacetime is dual to a quantum field 
theory on its conformal boundary. In this holographic 
setup, the radial dimension of space is emergent, but time 
is not: there is a global time which plays an entirely 
conventional role in the boundary theory. 

In this Letter, we use AdS/CFT to study a big crunch 
cosmology where the bulk singularity corresponds to the 
roll of a boundary scalar field to infinity in finite time. 
Several ingredients allow us to describe the dynamics in 
detail. First, the relevant coupling in the boundary the- 
ory is asymptotically free, rendering perturbation the- 
ory reliable. Second, the boundary evolution is smooth 
and ultralocal, even nonlinearly, as the singularity is ap- 
proached. Finally, the boundary has finite volume, so the 
homogeneous mode of the unstable field undergoes quan- 
tum spreading. With self-adjoint boundary conditions 
imposed at large field values, the wavefunction reflects 
off infinity with vanishing weight there, providing an au- 
tomatic UV cutoff on particle creation. Remarkably, the 
boundary theory develops nearly scale-invariant fluctua- 
tions as a consequence of its (slightly broken) conformal 
invariance, suggesting a new, natural mechanism for gen- 
crating cosmological perturbations. A detailed account 
of this work is given in Ref . [2| . 

Ten-dimensional type IIB supergravity compactified 
on S 5 can be consistently truncated, in several inequiv- 
alent ways, to five dimensional gravity coupled to a sin- 



gle scalar field tp, describing a quadrupole distortion of 
the S 5 .]$ For example, truncating to an S l O(5)-invariant 
scalar (p yields the Lagrangian 



1 



-R 



1 



1 



4 ^AdS 
vl/2 



(I5e 2w + 10e~ 4w - e 

with 7 = (2/15) 1 / 2 . The static AdSs vacuum occurs 
at ip = 0, where the scalar field mass saturates the 
Breitenlohner-Freedman (BF) bound m 2 = — 4i?^ 2 s . 
We are interested in time-dependent supergravity solu- 
tions which evolve into a big crunch singularity in the 
future (and possibly out of a big bang singularity in 
the past). In all solutions asymptotic to AdSs, with 
ds 2 



r 2 )dt 2 + dr 2 /(l + r 2 ) + r 2 dQ 2 }, tp tends 
, with a and (3 functions of the bound- 



to ar 2 In r + f3r 
ary coordinates, at large r. Specifying the dynamics re- 
quires a boundary condition on timclikc infinity relating 
a to [3. The supersymmetric, AdS-invariant choice is 
a = 0: but a — f/3, with / > 0, allows smooth, asymp- 
totically AdS data to evolve to a big crunch singularity 
in finite global time. 

For / = 0, the dual field theory is SU(N) gauge theory 
with Af = 4 supersymmetry on IR x S 3 . The bulk scalars 
saturating the BF bound couple to quadratic traces of 
Higgs scalars traceless on their £0(6) indices. For ex- 
ample, the S'0(5)-invariant bulk scalar tp, couples to 
O ~ A^Tr^ 2 ~ (l/5)Y?i= 2 where the trace is on 
SU(N) indices. The function a appearing in the bound- 
ary condition plays the role of a source for O in the 
boundary theory, while (3 corresponds to the expectation 
value of O. Boundary conditions with / > correspond 
to a deformation of the CFT: S -> S + J JO 2 . 1 i.e., a 
scalar potential driving O to infinity in finite time. Simi- 
lar cosmological duals were discussed in Refs. @ but the 
results were inconclusive because the dual 3d theories are 
poorly understood. [6j In the example chosen here, the 
dual 4d theory is a renormalizable theory with a pertur- 
bative regime. Furthermore, because of the "wrong" sign 
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of the deformation, the coupling / is asymptotically free. 
It follows that quantum corrections are under excellent 
perturbative control for large O, i.e., near the singularity. 
The large N effective potential reads Q 

V = -0 2 /\n(0/M 2 ) = -foO 2 . (2) 

Through dimensional transmutation, / has been replaced 
by an arbitrary mass scale M.[8j If fo is small at one 
value of O, it is even smaller for larger values, so one can 
be sure that quantum corrections do not turn the poten- 
tial around. At small 't Hooft coupling, we have control 
over the gauge theory. According to the AdS/CFT re- 
lation g t = §y M N = (RAds /^s) 4 , the bulk is then in a 
stringy regime. The supergravity regime corresponds to 
large 't Hooft coupling g t in the gauge theory. However, 
Witten's discussion [4[ indicates ([2]) is still valid. Hence 
the unbounded negative potential is a robust feature at 
small or large 't Hooft coupling. 

Since V is unbounded below, it is not obvious the field 
theory makes sense. A key point is that the S* 3 is finite, so 
the unstable homogeneous mode is quantum mechanical. 
Any initial wavepacket spreads instantaneously to infin- 
ity, so one simply cannot discuss the quantum evolution 
without imposing a boundary condition at large 0, cho- 
sen to ensure the Hamiltonian is self-adjoint. In quantum 
mechanics this construction is known as the "self-adjoint 
extension." [l(| We study its field theory implementa- 
tion and argue furthermore that perturbation theory is 
reliable for computing the amplitude to obtain almost all 
final values of the homogeneous mode. 

The self-adjoint extension is usually formulated [l(| 
as a boundary condition on energy eigenstates. We re- 
formulate it using the time-dependent semi-classical ex- 
pansion, which becomes exact near the singularity. [2| 
The Schrodinger equation is solved as an expansion in 
h, * = {A{x,t) + h...)e lS{x ^l h+ - where S(x,t) is the 
action for a complex classical solution ending at real x at 
time t. The initial condition for a Gaussian wavepacket 
centered on x c , with momentum p c , is x + 2i(Ax) 2 p/h = 
x c + 2i(Ax) 2 p c /h 1 where Ax is the initial spread in x. 
The imaginary part of the classical solution is thereby 
related to the quantum spreading of the wavefunction. 

For a potential like V(x) — — (l/4)A:c 4 , invariant un- 
der reflecting x through infinity in the complex x-plane, 
a wavefunction satisfying a self-adjoint boundary condi- 
tion may be constructed via the method of images. The 
"image" part of the wavefunction is obtained from a com- 
plex solution with (x c ,p c ) — > — (x c ,p c ) but the same final 
value of x. The complete semiclassical wavefunction is 

V a {x,t) = A 1 {x,t)e tSl / h + e- 2ta A 2 (x,t)e lS2/h , (3) 

where S± is the action for the first complex solution, S2 is 
that for the "image" solution, A4 are the corresponding 
amplitude factors and < a < ir is an arbitrary phase 
specifying the particular self-adjoint extension. At large 
x, we find the x-dependence to be 

ty a (x, t) —f A(t) x^ 1 cos + a) , (4) 



that of the self-adjoint energy eigenstates [T(|: thus ^f a 
lies in the domain of the self-adjoint extension. The phys- 
ical interpretation is that a right-moving wavepacket dis- 
appearing at infinity is always accompanied by a left- 
moving wavepacket appearing at infinity, with phase 
2a. To describe a bounce, we start with a Gaussian 
wavepacket rolling downhill, dominated by the first term 
in ([3]). The wavefunction reflects from infinity and rolls 
back uphill in a similar wavepacket, dominated by the 
second, "image" term. The phase a then nearly cancels 
in the final probability density. 

We can extend this technique to the dual field theory. 
We focus on the steepest negative direction of V, pa- 
rameterized as $1 = <f>U, with U a normalized SU(N) 
generator, Tr(U 2 ) = 1, and $ 4 = 0, i = 2 . . . 6. The 
effective potential becomes 

v ~ r aIs^- X ^^ X >P = V {N 2 ln(<f>/(NM))) = N- 2 l-\ 

( 5 ) 

ignoring numerical factors, with M = MN 1 ' 2 . We as- 
sume MRAds -C 1, so I is large for all <f> beyond the 
barrier. In this case, many quantities of interest can 
be expanded in inverse powers of I. Consider the clas- 
sical solution starting with zero kinetic and potential 

energy, at (j> ~ ^Ads- As rous to l ar g e val- 

ues, it approaches the scaling solution <f> ~ |t | 1 

plus l^ 1 corrections. Denning \ = d(j) 1 ^ 2 (/)~ 2 ~ 

A < ^ 1 ^ 2 _1 (1 + + ...), the solution is just \ = 1*1- 
Initially, x ~ RAdS- the time to bounce off infinity and 
return is ~ 2RAds- Now consider a Gaussian wavepacket 
in 0, centered on the above solution. Given an initial 
spread Afa, we can use Air^^ ~ R AdS A<fii ~ h/ A<f>i to es- 
timate the final spread, A\f ~ (Axf + R 2 AdS {Axi) 2 ) 1/2 ■ 

1 /2 

This is minimized by Axi, m in ~ A^ R A ds < Xc or 

A(f>i im i n ~ R-AdS ^ ^ c > wnere 0c is the center of the 
initial wavepacket and Xc the corresponding value of x- 
We shall write the spread of the initial wavepackets we 
consider as WA(j>i tm i n . 

To describe the bouncing wavepacket we must find 
the appropriate complex classical solution. For times 
t i=s +RAdSj the wavefunction is dominated by the image 
solution. Its initial condition is (f>+2ih~ 1 ir ( f > (A(f>) 2 = —<f> c , 
re-expressed in x- With the final condition X = r f> real, 
corresponding to <fi = (f>f, the solution is fully determined. 
To a first approximation I may be treated as constant 
and energy conservation leads to % = (1 + e^ 4 ) 1 / 2 with 
e a complex constant. Setting x = r + one finds r 
is close to the real, zero-energy scaling solution r = t, 
and s ~ (si/2) (l — (f/^/) 5 ) ■ Up to numerical factors, 
Si ~ W 2 (rf — rci), where rci is the classical value of r 
at the final time, and thus positive (negative) if the ar- 
gument of the wavefunction rj is ahead of (behind) the 
classical bouncing solution. 

Including the running of the coupling complicates 
the analysis because the logarithm introduces a branch 
cut at 4> = oo. However, imposing a UV cutoff A, the ef- 
fective potential V is meromorphic for \<f>\ > A and hence 
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the method of images still works. We carefully compute 
the relevant classical solution at large, finite A, following 
it onto the appropriate Riemann sheet before finally tak- 
ing A to infinity. The final wavefunction is A-independent 
and defines a completely unitary evolution. [2| 

In describing the bouncing homogeneous mode, we 
have so far ignored particle creation. We can check this 
is consistent by including the inhomogeneous modes in 
our two complex classical solutions. Each mode is placed 
in the incoming adiabatic vacuum state by imposing the 
initial condition that it is purely positive frequency. We 
then solve for the inhomogeneous modes to quadratic or- 
der, determining the amplitude of the outgoing negative 
frequency component. Finally, we compare the energy 
density in quantum-created particles with the kinetic or 
potential energy in the homogeneous mode. As we shall 
see, in a particular parameter regime, backreaction is 
negligible. 

When the modes of interest have wavelengths <C R-AdS, 
the curvature of the S 3 can be safely neglected. The 
equation of motion for fluctuations in <j> is 

84> = (^-k 2 + 6 X ~ 2 (1 + ^r 1 +...)j8<t>= -u 2 (t)5<t>. 

(6) 

The semiclassical expansion gives the wavefunction for 
the mode exactly to quadratic order, 

*(<5<M) ~ e i£ ?\ Sci = l(Rin/Rin){S4>) 2 , (7) 

with Ri n — ► e %UJt ,kt -C — 1, being the incoming positive 
frequency mode. The number of created particles at the 
final time tf is 

(-2iu(t f ))(R* n R in - Ri n it* n ) ' 

where w(tf) is real. Now, a very important fact comes 
into play. Near the singularity, the complex solution for 
the homogeneous mode is well-approximated by x = t + 
ie, where e = s,-/2 is a constant of order W 2 \ 1 J 2 RAdSi IOT 
typical final values of the homogeneous background field 
(j>f. For most values, therefore, the dangerous x~ 2 term 
in © is actually bounded and smooth at t = 0. Hence 
quantum particle production is exponentially suppressed 
for k > lei" 1 ~ W- 2 \- 1/2 R A l ds ~ W^Nl^R^. 

At leading (zeroth) order in the solution Ri n of ([6]) 
is (t + ie) 1 ^ 2 H 5 / 2 (k(t + ie)) with the property that a pos- 
itive frequency incoming mode evolves to a positive fre- 
quency outgoing mode. This behavior was first noticed in 
field theory on compactified Milne spacetime, where e was 
used to define a choice of analytic continuation. Here, 
e is a physical parameter determined by the final value 
of the homogeneous "background," i.e., the argument of 
the wavefunction, and except when e vanishes, the evolu- 
tion of linearized fluctuations is unambiguous. Solving to 
next order in Z _1 , we find the positive frequency incoming 
mode evolves to ~ e lkt — ine~ lkt e~ 2k \ e \ / ln(fc/M) in the 



outgoing adiabatic regime, kt ^> 1, t ~ RAdS *C M , 
giving the energy density in quantum-created 5cj> particles 

f d 3 k it 2 _ 4fc | £ | A^ 4 
PcM - J (2^ ln 2 (fc/M) e ~ W8R\ ds ' (9) 

where the logarithm enters as the value of I when a mode 
leaves the adiabatic regime, when \kt\ ~ 1. 

We now consider quantum production of other modes. 
Assuming $! points in general position in SU(N) space, 
i.e., SU(N) is broken to [/(l)^ -1 , N gauge, Higgs and 
Fermi particles remain light. The other Higgs and gauge 
particles acquire masses M ~ gy m4* with M/M 2 ~ 
X^/ffYM ~ {gtlN)- 1 ' 2 < 1 at large I and N, so 
the corresponding positive frequency mode evolution is 
adiabatic and particle production is exponentially sup- 
pressed. 

The light Higgs particles acquire mass from the term 
A Tr($?)Tr($f), i = 2,..., 6. Their mode functions 

obey = -fc 2 $« - (2/5)(l + r 1 + 0(l- 2 ))^/ X 2 - To 
zeroth order in I , Ri n is again a Hankel function and 
yields no particle production. To first order in l^ 1 , up 
to a numerical factor we find an energy density N times 
greater than 

The massless gauge bosons have no classical cou- 
pling to 4>. However, a coupling arises at one loop. 
The effective Lagrangian for the massless gauge fields 

is -a/^)F 2 (l + g 2 YM E l =s,F,v^HM 2 / f i 2 )), with 
fj, a renormalization scale and the trace taken over the 
heavy scalars, fermions and gauge bosons, with mass 
matrices .Mi.[l2| In the undeformed theory, the correc- 
tion vanishes by supersymmetry. The deformation al- 
ters M% ~ g 2 <& 2 —> g 2< &\ + A<^$ 2 , producing a finite, 
/i-independent correction ~ N~ 1 l~ 1 F 2 . The equation 
of motion for the massless gauge boson modes yields a 
mode mixing coefficient ~ N~ 1 e~ 2k ^ / ln(fc/M) 2 across 
the bounce. Compared to ©, the produced energy den- 
sity is suppressed by one power of N and two powers of 
\n{k/M). 

A conservative criterion for backreaction to be small 
is that the homogeneous mode should keep enough 
energy to roll back up the hill to nearly its start- 
ing value, (j) s tart ~ R Ads N Q- n (<f>start/NM)) 1/2 w 
R~AdsN\ ^(MRAds)] 1 ^ 2 ■ The minimal value of <j) reached 
after the bounce is obtained by setting V(<p m i n ) = —p c - 
Demanding that 4> min - 4> s tart < 4>start in order to 
have small backreaction, one finds the condition p c <C 
A^ 2 | \n{M RAds)\/ R\ds'i as tne l &r g es t contribution to p c 
comes from light Higgs particles with p c ~ N 5 /W 8 R Ads , 
this becomes 

\]n{MR AdS )\>N 3 /W 8 . (10) 

So, for example, if we choose W ~ N 3 ^ 8 , backreaction 
is under control as long as \\n(MRAds)\ ^S> 1, for fi- 
nal values of <fi carrying most of the probability. This 
value of W corresponds to a spread in the time delay 
e ~ A^ 1 / 4 ! ln(MRAds)\~ 1 ^ 2 RAds, parametrically small 
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compared to the duration ~ 2Ra<is of the bounce. Hence 
our wavepacket is well-described by a classical bouncing 
solution, with small backreaction. The fraction of the 
probability associated with final cj> values for which the 
backreaction is large (and hence our calculation breaks 
down) is ~ \ln(MRAds)\~ 1 ' 4 '- Hence we conclude that 
when the logarithm is large, the most probable outcome 
is a bounce. 

While the main emphasis of our work is to establish 
the possibility of a nonsingular bounce, the fluctuations 
in the rolling field <f> have a further remarkable conse- 
quence: they give rise to a nearly scale-invariant spec- 
trum of stress-energy perturbations in the dual theory. 
As <j> rolls, its fluctuation S(f> acquires an increasingly 
tachyonic mass, causing modes of successively higher k 
to cease oscillating and enter a "growing mode" solution, 
just as in the ekpyrotic mechanism. [h| In that case, the 
scale-invariant spectrum generated is the result of the 
classical scale-invariance of a scalar field with a nega- 
tive exponential potential under rescaling spacetime co- 
ordinates and shifting the field 0-^0 + const. Here, 
the field theory Lagrangian is also unstable and classi- 
cally scale-invariant, with the scaling symmetry rescal- 
ing <f> according to its canonical dimension. The classi- 
cal instability leads naturally to a nearly scale-invariant 
spectrum of stress-energy fluctuations on the boundary. 
Quantum effects break the scale symmetry but with small 
anomalous dimensions because is asymptotically free. 
A straightforward but lengthy computation of the two- 
point correlator for the "improved" stress-energy tensor 
yields for the dimensionless fluctuation 5 = Sp/(P + p), 
with p = T 00 and P = T il /3 as usual, the two-point cor- 
relator (5(r,t)S(0,t)) ~ N- 2 {\n Mr)- 2 (In Mt^Fir/t) 
with F a dimensionless function. The only violations of 
naive (engineering dimensions) scale-invariance are due 
to the logarithms and they are small at weak coupling. 
The stress-energy perturbations are (i) naturally small, 
being suppressed by powers of 1/N and A^, (ii) approx- 
imately scale-invariant, with a slightly red tilt due to 
asymptotic freedom, (iii) scalar, (iv) adiabatic and (v) 



nearly Gaussian in character, to leading order in A^ and 
1/N. This is a tantalizing result. To be sure, we have 
only so far computed boundary correlators and hence not 
the quantities of direct interest in the bulk cosmology. 
However, since the only scale which enters the translation 
is R-AdS and we are interested in modes with k ^> R^dg, 
and since scale-invariance is "holographic" in the sense 
that the fluctuations on a planar subspace of a space 
with scale-invariant fluctuations are scale- invariant, it is 
plausible that the bulk cosmological perturbations will 
turn out to be scale- invariant. (For related calculations 
in 0(1, 3) -invariant backgrounds see Ref. [3.) 

A more realistic cosmology requires a setup with a 
four dimensional bulk, like AdSi x S 17 .^ The conformal 
boundary is 1R x S 2 and the dual theory is believed to 
be the IR limit of a 2 + 1-dimensional, large N gauge 
theory. As in our case, the dual theory is expected to be 
classically conformal-invariant and unstable: the model 
dynamics proposed in Ref. Q consists of a scalar field <fi 
with a deformation potential —f<fi 6 . In this model, / is 
again asymptotically free. While precise calculations are 
currently beyond reach, it is not unreasonable to expect 
qualitatively similar behavior to that found here. If so, 
and if the model can be embedded in a bigger picture al- 
lowing for a positive dark energy phase, this mechanism 
may provide a rather compelling explanation for the ori- 
gin of the observed density variations in the universe. 
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